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We measure the elastic and dissipative responses of a multicontact interface, formed 
between the rough surfaces of two contacting macroscopic solids, submitted to a 
biased oscillating shear force. We evidence that beyond a linear viscoelastic regime, 
observed at low shear amplitude, the interface response exhibits a dissipative com- 
ponent which corresponds to the onset of frictional dissipation. The latter regime 
exists whereas the tangential force applied, far from the nominal static threshold, 
does not provoke any sliding. This result, akin to that of Mindlin for a single con- 
tact, leads us to extend his model of 'microslip' to the case of an interface composed 
of multiple microcontacts. While describing satisfactorily the elastic response, the 
model fails to account quantitatively for the observed energy dissipation, which, 
we believe, results from the fact that the key assumption of local Coulomb friction 
in Mindlin's model is not legitimate at the sub-micrometer scale of the microslip 
zones within microcontacts between surface asperities. 
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1. Introduction 

The frictional response of the contact between two macroscopic solids submitted to 
a shear force is commonly described, in the framework of Amontons-Coulomb's law, 
in terms of: (i) a static force threshold, F s , below which no relative displacement is 
supposed to occur, and (ii) a dynamic friction coefficient defined when stationary 
sliding is established. 

However, it is known that frictional dissipation in mechanical contacts starts to 
build up for shear forces lower than the nominal static threshold. This behaviour, 
which is important in mechanical engineering, for instance in problems of fret- 
ting or damping in structural joints (see e.g. Goodman 1959, Olofsson 1995), also 
presents a fundamental interest related to the understanding of the microscopic 
mechanisms responsible for macroscopic friction. This issue was first extensively 
studied by Mindlin et al. (1951), Johnson (1955), Courtney-Pratt & Eisner (1956) 
and Goodman & Brown (1962), who evidenced that the displacement response of 
Hertzian macroscopic contacts submitted to an oscillating tangential force of am- 
plitude F -C F s exhibited a hysteresis loop attributable to an interfacial dissipative 
process. In order to account for this energy loss, Mindlin et al. (1951) proposed 

f Corresponding author, e-mail: lionel.bureau@college-de-france.fr. Present address: LPMC, 
College de France, 11 place Marcelin Berthelot, 75231 Paris cedex 05. 



Article submitted to Royal Society 



TgX Paper 



2 



L. Bureau, C. Caroli and T. Baumberger 



a model of 'microslip' within the contact zone, based on the following description 
(Cattaneo 1938): wherever in the contact zone the tangential (r) and normal (a) 
local stresses obey t < jia (with fi the friction coefficient) no relative tangential 
displacement occurs. Otherwise, shear provokes slip so that, in the corresponding re- 
gion, the stresses satisfy locally a Coulomb friction law: r = pta. It is then predicted 
that the sheared contact is composed of a circular non-sliding zone surrounded by 
an annular slipping region whose width increases with the applied shear force. 

Still, Mindlin's description of incipient sliding raises a question: down to which 
length scale is such an assumption of local friction valid ? Indeed, macroscopic solids 
generally exhibit rough surfaces which, when brought into contact, form a multi- 
contact interface (MCI) , i. e. an interface composed of a dilute set of microcontacts 
between asperities. When addressing the problem of friction at such an interface, 
one may therefore wonder whether the use of a local friction law is legitimate when 
dealing with the micrometer-sized contacts between surface asperities. This in turn 
raises the question of the spatial scale of the elementary dissipative events respon- 
sible for solid friction, a cutoff length below which a local friction law cannot be 
meaningful. 

Recently, experiments performed on such MCIs showed that (Berthoud & Baum- 
berger 1998): 

• for tangential forces F <C F s , the pinned interface responds elastically, via 
the reversible deformation of the load-bearing asperities. In this regime, the 
shear stiffness k of the interface can be measured and is well accounted for 
within the framework of Greenwood's model of contact between rough surfaces 
(Greenwood & Williamson 1966). 

• for F < F s , a creeplike irreversible sliding of the solids occurs. 

A detailed study of this regime of incipient sliding, below the nominal static 
threshold, should therefore provide information about the physical processes un- 
derlying frictional dissipation. In order to achieve the force control and the dis- 
placement resolution required to perform such a study, we have developed an ex- 
perimental setup which allows to probe the response of a multicontact interface to a 
biased oscillating shear force, of amplitude F ac about a finite mean value F^c, while 
the maximum force applied Fd c + F ac = F max < F s (Baumberger et al. 1998). We 
thus showed that the displacement response of a macroscopic slider submitted to 
such a harmonic tangential force exhibited three different regimes, depending on 
the amplitude F ac (Bureau et al. 2001). These are illustrated on figure 1, which 
shows the response to a shear force modulation of slowly increasing amplitude: in 
region (i), at small F ac , the center of mass of the slider oscillates about a con- 
stant average position, which means that no irreversible sliding occurs; in region 
(ii), corresponding to higher shear amplitudes, the slider enters a creeplike regime 
where the slipped distance increases continuously, up to a final regime of abruptly 
accelerating motion (region (iii)). 

In the present article, we report an experimental study which focuses on the first 
regime of small shear amplitudes, in which the interface is pinned, i.e. oscillates 
about a fixed position. In section 2, we briefly present the principle of operation of 
the experimental setup. By measuring the in-phase and out-of-phase components 
of the displacement response of a MCI submitted to a harmonic shear force, we can 
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Figure 1. Time-plot of the instantaneous (line) and averaged (•) displacement response 
of the slider to a biased oscillating shear force of increasing amplitude. On the right scale 
is plotted the reduced shear amplitude (o) 7 ac (£) = F ac (t)/W , where W is the normal 
load. The biased used for this experiment is 7d c = Fdc/W = 0.36 and the frequency of 
oscillation / = 80 Hz. 

probe accurately both the elastic and dissipative responses of the interface (section 
3). We evidence that: 

• the shear stiffness evolves with the 'age' of the interface, i.e. with the time 
elapsed since the solids were brought into contact, in agreement with the 
creep ageing of the load-bearing asperities already identified for such MCIs, 
which results in the slow logarithmic increase of microcontact radii (Dieterich 
& Kilgore 1994, Berthoud et al. 1999a), 

• frictional dissipation appears at shear forces well below the nominal static 
threshold, while no gross sliding is detected yet. 

This latter point leads us to propose an extension of Mindlin's model to the case 
of a multicontact interface, within the framework of Greenwood's description for 
the contact of rough surfaces. We find that such a model of microslip within the 
contacts between asperities largely overestimates the energy dissipation, which we 
believe points out that a local description of friction is no longer valid at the sub- 
micrometric scale of the slip zones within the microcontacts forming multicontact 
interfaces. 

2. Experiments 

(a) Samples 

The multicontact interface studied here is formed between a slider and a track of 
commercial grade poly(methylmethacrylate) (PMMA). PMMA is a glassy polymer 
at room temperature (T g ~ 120°C), of dynamic shear modulus G' = 2 GPa and 
loss angle 5 such that tan 5 = 0.1 at 100 Hz (Ferry 1980), and of quasi-static Young 
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modulus E ~ 3 GPa, Poisson ratio v = 0.44 and hardness H ~ 300 MPa (Berthoud 
et al. 1999b). 

The nominally flat surfaces of the slider (20x20 mm 2 ) and of the track (25x30 
mm 2 ) are lapped with an abrasive aqueous suspension of SiC powder (mean particle 
size 23 /nm). This leads to a rms roughness of the samples R q = 1.3 /im, as previously 
characterized (Berthoud & Baumberger 1998). 

Using Greenwood's result (Greenwood & Williamson 1966), assuming that: 
(i) the summit heights of asperities follow an exponential distribution of width 
s = 1.3 /im, and (ii) their mean radius of curvature (3 = 20 /im (taken, as a con- 
servative value, of the order of magnitude of the abrasive particle size), one can 
estimatef, under the normal load W ~ 2 N, the number N of microcontacts: 
N ~ W/(^/^j3E*s) ~ 100 (where = £/[2(l - ^ 2 )]), and their mean radius 
a ~ \fsj3 — 5 /im. 

(b) Experimental setup: an inertial tribometer 

The experimental setup, extensively described elsewhere (Baumberger et al. 
1998) is based on the following principle. The track, on which the slider sits under 
its own weight mg, is first inclined at a given angle 6 such that the ratio of the 
tangential (Fd c ) to normal (W) load Fd c /W = tan# <C fJ, s - A harmonic motion, 
of controlled amplitude and frequency, is then imposed to the track, which results 
in an inertial oscillating shear force acting on the center of mass of the slider (see 
figure 2). In order to avoid torques and tilt motion, the slider sample is clamped in 
a metal part specially designed such that the center of mass of the slider is located 
in the plane of the interface. 

The harmonic shear force on the interface thus reads, in the low frequency limit 
where the slider responds quasi-statically (see Baumberger et al. 1998), F(t) = 
F ac cos(ujt) with F ac = mT, where T is the imposed acceleration amplitude of the 
track. We define jdc — Fd c /W and j ac — F ac /W. These control parameters can be 
set in the ranges < jdc < 0.58 and < j ac < 0.6. The frequency of the oscillating 
tangential force can be chosen between 15 and 100 Hz, this upper limit ensuring 
the quasi-static condition for the slider motion. 

We measure, in response to this excitation, the displacement x(t) of the slider 
relatively to the track, by means of a capacitive gauge. Its signal is sent to a lock-in 
amplifier, which allows to detect, within 1 nm, the in-phasc (x ) and out-of-phasc 
(xqo) components of the displacement with respect to the harmonic input. We thus 
have access to the elastic and dissipative responses of the multicontact interface. 

(c) Reproducibility 

The scattering of the results thus obtained depends crucially on the way the 
interface is prepared. We have tested two different protocols: 

f This estimation assumes that contacting asperities arc deformed elastically. Actually, Green- 
wood's plasticity index ip = (E/H)y^s/f3 ~ 2.5 for our surfaces, which indicates that most con- 
tacting asperities have started to yield plastically. However, since rp is still of order unity, the set of 
asperities is in fact in an elastic-plastic state of deformation (Berthoud et al. 1999b), far from the 
fully plastic limit. We therefore use, for our rough estimate of N, the expression corresponding to 
the elastic limit (evaluating in the fully plastic limit would lead to a number of microcontacts 
ip times larger). 
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Figure 2. Principle of operation of the inertial tribometer: the slider on the inclined track 
is submitted to a constant tangential force Fa c = nig sin 8 on which is superimposed a 
harmonic shear force of amplitude mT, where F is the imposed acceleration amplitude of 
the track. We measure the displacement x of the slider with respect to the track. 

(i) after the slider is put into contact with the inclined track, a time t w is waited 
during which the interface is left under constant normal and tangential loads. At 
the end of this waiting time, j ac is turned on for a time <C t w during which we 
measure Xq and Xgo ■ When performing several such experiments in the same nominal 
conditions, the relative dispersion observed is on the order of 25%. 

(ii) the second protocol consists, as soon as the interface has been created, in 
applying to the slider a large amplitude harmonic shear force in order to make it 
slide a few micrometers in the direction of Fdc- The oscillating shear force is then 
suddenly stopped and a time t w is waited, after which measurements are performed. 
The dispersion observed when using this second protocol is reduced to 11%. We 
will therefore present, in the rest of the paper, results that have been obtained on 
interfaces prepared this way. 

This effect of the interface preparation on the reproducibility is illustrated on 
figure 3. 
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Figure 3. Influence of the protocol of interface preparation on the reproducibility. In-phase 
response a;o(7ac) obtained with protocol (i) (o), and with protocol (ii) (•). Four sets of 
measurements are plotted, performed at / = 60 Hz, with 7<j c = 0.27 and after t w — 300 s. 
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To understand, at least qualitatively, the origin of this effect, note the following. 
When the slider is put into contact with the inclined track, the only mechanical 
condition which constrains the state of the MCI is that the total tangential force 
F = WtsmO. This is obviously insufficient to define uniquely the distribution of 
shear forces on the various microcontacts. Hence, a huge number of local configu- 
rations are possible: the interface is a highly multistable system. On the contrary, 
sliding produces a reproducible distribution which is maintained during the elastic 
recoil following a stop (Caroli & Nozieres 1996, de Gennes 1997). Preparing an 
interface by interruption of sliding can therefore be expected to improve the degree 
of reproducibility. The remaining scattering of 11% is indeed in agreement with the 
statistical dispersion ~ N~ x / 2 ~ 10% due to the finite number of microcontacts 
estimated above. 



3. Results 

We present on figures 4 and 5 typical results obtained for the in-phase (x ) and 
out-of-phase (xgo) components as a function of the reduced shear force amplitude 
7 ac . We do not observe any dependence of the elastic or dissipative response on the 
level of average tangential force jdc («-e. on the angle of inclination of the track, up 
to -fdc — 0.36), and do not note any frequency dependence in the explored range 
15-100 Hz. 
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Figure 4. In-phase displacement amplitude xo as a function of 7 ac , at / = 80 Hz and 
after t w = 600 s. (•): 7dc = 0, (o): 7<j c = 0.09, (□): -y dc = 0.27. 



Both xo(7 ac ) and a;go(7 ac ) exhibit a linear regime at low shear amplitude, up 
to 7ac — 0.1 above which the interface responds non-linearly. For all the results 
presented in this article, in the range of j ac explored, no sliding is detected within 
the experimental resolution < 20 nm. 
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Figure 5. Out-of-phase displacement amplitude xgo as a function of 7 ac , at / = 80 Hz 
and after t w = 600 s. (•): j dc = 0, (o): 7dc = 0.09, (□): -y dc = 0.27 



In the linear regime, we measure the interfacial shear stiffness, which reads 
k = IucW/xq. Moreover, we know that, for a multicontact interface, this stiffness 
varies proportionally to the normal load: k — W/X, where A is an 'elastic length' 
which lies in the micrometer range (Berthoud & Baumberger 1998, see also section 
4a for more details). The slope of xo(7ac) is thus A = dxo/d7 ac . We measure, after 
a waiting time t w = 300 s, A = 0.26 ± 0.015 /an. 

This length can also be determined by analyzing the frequency response of the 
system. Indeed, we expect that the response of the slider, of mass to, sitting on 
the elastic foundation of stiffness k formed by the set of load-bearing asperities, 
exhibits a resonance at a circular frequency u a — yjn/m = y/g/X (with W = mg), 
as already observed by Sherif & Kossa (1991). When looking at the spectral response 
shown on figure 6, one clearly identifies a peak at f = 1000 Hz which corresponds 
to that resonance and leads to a value of the elastic length A = 0.25 /tm, in excellent 
agreement with the value reported above. 

Let us now try to identify the origin of the dissipation in this regime of very small 
shear amplitudes. We note that the inverse of the quality factor 1/Q = 0.03 of the 
resonance of figure 6 is comparable with the tangent of the loss angle tan 5 ~ 0.05 of 
bulk PMMA at 1kHz and T = 300 K (Ferry 1980). Besides, the ratio x 90 /x = 0.18 
is constant and on the order of tan S w 0.1 at / < 100 Hz. This leads us to 
attribute the observed dissipation to the viscoelastic losses within the bodies of 
contacting asperities - more precisely, for each of them, within the volume of order 
roughly a 3 (with a the mean contact radius) in which stresses concentrate. The 
discrepancy observed between xgo I and the loss angle measured on bulk samples 
can be assigned to the fact that these micrometric volumes lie in the interfacial 
region, and hence most probably present mechanical properties slightly different 
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Figure 6. Spectral response: the slider is submitted to an acceleration of constant amplitude 
and variable frequency. The ratio of the displacement amplitude Aa; to the acceleration 
amplitude is plotted, in arbitrary units, as a function of /. The bias 7 dc = 0. 

from those of the bulk. Indeed, our method for surface abrasion (see §2a) makes 
use of water, known as a plasticizer of PMMA. 

(b) Ageing 

When performing measurements in the linear regime at various waiting times 
t w , we note that the interfacial shear stiffness, or equivalently the elastic length 
A, evolves slowly with t w : the longer t w , the lower A, i.e. the higher the stiffness. 
This ageing of the interface is illustrated on figure 7. The elastic length decreases 
quasi-logarithmically with the age of the interface: 

X(t w ) = A -/?aMW*o) (3.1) 

with the reference time t — 1 s, A = 0.33 and the logarithmic slope (5\ = 1.07 x 
TO" 2 . 

(c) Non-linear regime: interfacial dissipation 

At shear amplitudes -f ac > 0.1, the in-phase and out-of-phase components of 
the response increase non-linearly. In this regime, while the bulk response of the 
asperities remains linear (one can estimate a mean shear strain as the ratio of x to 
the mean contact size a, which stays lower than 2%), the ratio x 90 / x o 1S n °t constant 
anymore, hence the energy loss cannot be attributed to bulk viscoelasticity only. 

We show on figure 8 the evolution, with shear amplitude, of the non-viscoelastic 
part of the dissipative response, which we define as xgo — xo tan 6, using for tan 5 
the value xgo/xo = 0.18 determined in the linear regime. 
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Figure 7. Decrease of the elastic length A with the waiting time t w . f = 80 Hz, 

7dc = 0.18, j ac = 0.04. 
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Figure 8. Non-viscoelastic component of the dissipative response, £90 — :rotan<5, as a 
function of "/ ac - f = 80 Hz and t w = 600 s. (•): 7<j c = 0, (o): j dc = 0.09, (□): 7<j c = 0.18, 
(x): ldc = 0.27, (o): 7dc = 0.36. 



4. Discussion 

(a) Interfacial shear stiffness 

Let us consider the case of an infinitesimal shear modulation. In this limit, 
the shear stiffness of a single contact is simply given by the Mindlin expression 
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calculated in the absence of slip, namely Qa, with a the contact radius and Q = 
AG/ (2 — v), where v and G are respectively the Poisson ratio and shear modulus 
of the contacting materials (Mindlin 1949). The shear stiffness of the multicontact 
interface then reads (Berthoud & Baumberger 1998): 



with N the number of microcontacts and a their mean radius (this result is rederived 
in the Appendix). An important feature of Greenwood's description of the contact 
between rough surfaces is that the number of microcontacts varies linearly with the 
normal load W, whereas their mean size is independent of W. So, the interfacial 
stiffness is proportional to the load, i.e. n — W/X, with an elastic length reading: 



Using for the PMMA shear modulus its quasi-static value G ~ 1 GPa, Berthoud & 
Baumberger found from equation 4.2 an elastic length A ~ 1 /zm, in agreement with 
their quasi-static measurements of the stiffness during loading- unloading cycles. 

Our results, however, lead to values of A ~ 0.25 fim much smaller than those 
previously reported. This marked difference may have two distinct origins: 

• First, we clearly see from our experiments that the linear regime of inter- 
facial response corresponds to elastic displacements of the slider of at most 
20 nm (see figure 4). Such a resolution could not be achieved in the previ- 
ous quasi-static experiments, and the stiffness measured in that study was 
most probably underestimated, due to non-linear effects, which thus led to 
overestimated values of A. 

• Moreover, note from expression 4.2 that the elastic length is inversely pro- 
portional to the shear modulus. We thus expect the elastic response of the 
multicontact interface to be governed, in our experiments, by the dynamic 
modulus at the excitation frequency, i.e. G' ~ 2 GPa (see §2a). 

With this latter value for the dynamic modulus, along with those for the hardness 
and the surface characteristics given in section 2a, we obtain A ~ 0.39 /im. The elas- 
tic length that we estimate from the elastic properties of PMMA at the excitation 
frequency is therefore in good agreement with our experimental results. 



We now address the question of the time-dependence of the elastic length. It is 
well established, since the pioneer work of Bowden & Tabor (1950), that the real 
area of contact (S r ) between rough surfaces is in general a very small fraction of 
the nominal area. The normal stresses on the load-bearing asperities are thus on 
the order of the yield stress of the contacting materials, which results in bulk plastic 
creep of these asperities. As a consequence, the real area of contact slowly increases 
with the 'age' of the interface, i.e. with the time of contact between asperities, as 
unambiguously evidenced by Dieterich & Kilgore (1994). From an extensive study 
of the mechanical properties of polymer glasses, Berthoud et al. (1999a) have shown 



k ^ NQa 



(4.1) 




(4.2) 



(b) Ageing 
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that creep of the load-bearing asperities results in a quasi-logarithmic growth of X r : 



I f mln(l + ^ 



(4.3) 



with m and i c two material parameters that can be identified from bulk mechanical 
tests. For PMMA at room temperature, m e [0.04 — 0.05], and a higher boundf for 
the cut-off time is t c < 5 x 10~ 3 s. 

From expression 4.1 of the shear stiffness k = NQa = W/X, and with the real 
area of contact S r = No 2 , we find that the elastic length reads: 



A = 



W 



G^NT, r0 [l + mln(l + t w /t c ] 
To first order in m, this expression reads: 



A = 



W 



Ao 



™ln(l + ^ 

2 V to 



(4.4) 



(4.5) 



We thus expect the elastic length to decrease quasi-logarithmically with t w , 
with a logarithmic slope of m/2. Indeed, when fitting the data of figure 7 with 
an expression of the form A = Ao [1 + £hi(l + t w /t c )], leaving £ and t c as free 
parameters, the best fit is obtained for £ = 0.024 and t c = 10~ 3 s, these values 
being in full agreement with the expected ones. The observed dependence of the 
interfacial shear stiffness on the waiting time t w thus results from the creep ageing 
of the microcontacts. 

It is interesting to note that up to now, this mechanism of interfacial ageing 
has always been characterized through the time-dependence of the static friction 
threshold (see Berthoud et al. 1999a, and references therein), which is a 'destructive' 
method in the sense that the set of load-bearing asperities is renewed by sliding 
when the measurement is performed. On the contrary, our low-amplitude oscillating 
shear experiments provide a way to probe accurately the slowly evolving viscoelastic 
response of a given set of microcontacts, without any macroscopic sliding at the 
interface. 



(c) Non-linear elasticity and energy dissipation: extension of Mindlin's model 

Figures 4 and 5 clearly show that, for j ac > 0.1, the interfacial elastic response 
becomes non-linear, while no gross sliding is observable. In the same regime, a 
non-linear dissipative response develops on top of the linear term attributable to 
bulk viscoelasticity (figure 8). This contribution must therefore be considered as 
resulting from interfacial dissipation proper. 

The decrease of the 'local' interfacial stiffness, dj ac /dxo, with increasing shear 
amplitude may be interpreted qualitatively as follows. The diameters of the mi- 
crocontacts which form the interface are statistically distributed about the average 
value a. A finite shear necessarily leads to destroying the smaller ones. The larger 

f The cut-off time t c coufd not be determined accurateiy at room temperature, and was inferred 
from the velocity dependence of the friction force (Baumberger et al. 1999, Bureau 2002). 
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the shear amplitude, the smaller the number of microcontacts which are still able 
to sustain the stress, hence a decreasing stiffness. 

In order to describe quantitatively this regime, we now extend Mindlin's de- 
scription to the case of a multicontact interface as follows: 

(i) The contact between the two rough surfaces is described a la Greenwood, 
with the assumption of an exponential distribution of summit heights, and of elastic 
deformation of asperities. We believe the latter assumption to be inessential: indeed, 
we saw in §2a that though contacting asperities are in an elastic-plastic state, the 
value of the plasticity index ip is of order unity, which suggests that the normal 
stress profile in microcontacts is still close to the Hertz profile (Johnson 1985). 

(ii) Mindlin's results give, for a given microcontact, the expression of the tan- 
gential force associated with a remote shear displacement x. This displacement, 
equal to that of the center of mass of the slider, is common to all microcontacts. 

(iii) For any finite x, there always exists a set of small microcontacts which 
are completely slipping. For these, the tangential force is saturated at its constant 
maximum value nw, where w is the normal load on the microcontact. 

(iv) The tangential force on the slider is simply the sum of those on the various 
microcontacts. 

The detailed calculation is performed in the Appendix. It yields the following 
results: 



x n = 2/zA 



'Jac . i 'Jac \ | / 7ac \ //-if ^ ac 



■ . o ■ 



(4.6) 



(4.7) 

with A the elastic length defined above. The local friction coefficient of Mindlin's 
model, /i, is our single fitting parameter. 

Figure 9 shows the best fit thus obtained for the elastic part of the response 
x , which is seen to be excellent. It corresponds to fj, = 0.49. On the other hand, 
from the response to a linear ramp of shear amplitude (see figure 1), we have 
estimated the ('global') static friction coefficient /j, s as corresponding to a threshold 
of accelerated sliding. We thus find /j, s — 0.59 ±0.03 (Bureau et al. 2001). This can 
only be considered as a rough estimate, in view of the arbitrariness in the definition 
of the threshold. The agreement therefore appears quite satisfactory. 

At this stage, we have determined all the parameters of the model, and are thus 
in a position to truly check its validity by comparing its prediction for xgo(7oc) 
with the experimental results. The dissipation calculated from equation (4.7) with 
/i = 0.49 is plotted on figure 10: it is seen to be much larger than the experimental 
one. In order to get a decent fit of the data, we have to use a value of /x as large 
as 1.3, i.e. much larger than any value of the friction coefficient ever reported for 
PMMA. The question therefore arises of what might be the physical reason for such 
a discrepancy which affects primarily the dissipative part of the response. 

In Mindlin's model, dissipation results from the slip at the periphery of the 
contact. The inner radius of this annulus is, for an average microcontact, c = 



^90 



4/xA 

7T 



2 ^lnfl 7< " 



2^A 
"3V 
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Figure 9. (•): experimental data for xo(7 ac ) at / = 80 Hz and 7 dc = 0.09. ( — -): extended 
Mindlin's model with A = 0.23 /im and u — 0.49. 



0.01 



to 

§ 0.004 




-0.002 



Figure 10. (•): experimental data for a?9o(7oc) at / = 80 Hz and ~fdc = 0.09. ( — ): extended 

Mindlin's model with A = 0.23 /tm and u = 0.49. ( ): extended Mindlin's model with 

A = 0.23 /im and /i = 1.3. 



1/3 

a [1 — f/(fj,w)] . For a = 5 /mi and taking f/w = 0.2+ we find, for \x = 0.49 that 
the width of the annulus a — c = 750 nm. Over this distance, the shear stress varies 
from its maximum value to zero at the edge of the contact. We have thus tacitly 

f We use for the ratio f/w a typical value of the macroscopic ratio F/W. This amounts to 
assuming that JV identical microcontacts of size a bear the same fraction of normal (w = W/N) 
and tangential (/ = F/N) load. 
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assumed that the Coulomb law is valid on a spatial scale much smaller than this 
width. 

It is now well documented that frictional dissipation results from the dcpinning 
of structural elements located within the adhesive junctions between load-bearing 
asperities. The typical size of these elements is found to be nanometric (Nakatani 
2001, Baumberger et al. 1999, Bureau et al. 2002). The friction force, as it is usually 
defined, is an average over the dynamics of a large ensemble of such elements. So, a 
reasonably meaningful friction coefficient cannot be defined on a scale smaller than, 
say, a hundred nanometers. This must be understood as a cut-off length below which 
the Mindlin stress profile probably becomes inaccurate. The above estimate of a — c 
therefore suggests that multicontact interfaces with micrometric roughness might 
be out of the range of quantitative applicability of Mindlin's model. 



Appendix A. Extension of Mindlin's model 

(a) Mindlin's results 

We first recall the results derived by Mindlin for the contact of identical elastic 
spheres submitted to a tangential force (Mindlin et al. 1953, Johnson 1985). When 
the contact is first loaded from zero, the relationship between the remote tangential 
displacement x and the applied shear force / reads: 



3(2 — v)fj,w 
8Ga 



1 



1-1- 

fJ,W 



2/3' 



(Al) 



with G the shear modulus, v the Poisson ratio, a the Hertz radius of contact, w the 
normal load and \i the local friction coefficient. 

If the shear force is then decreased, after having reached a maximum value / mol , 
the displacement x\^ in this unloading phase is: 



3(2 — v)nw 
8Ga 



2 1- 



fit 



1 



2/iw 



2/3 



1 - 



In 



2/3 



1 



By symmetry, if the shear force is then reversed from a value — /„ 
placement x y(f) = —x\(—f). 



(A 2) 
c the dis- 



(b) Contact between rough surfaces 

The contact geometry is that of two rough surfaces of identical rms roughness a. 
We shall consider, within the framework of Greenwood's model, the case of contact 
between a rigid, ideally smooth, reference plane and a composite rough surface 
whose summit heights are distributed exponentially: (j)(z) = s~ 1 exp(—z/s), with 
s = \f2a (Bcrthoud & Baumberger 1998). The coordinate z is normal to the mean 
plane of the random surface. 

The effective elastic modulus of the dcformable material is defined as E* = 
E/[2{1 — v 2 )], and the equivalent radius of curvature at the tip of asperities for the 
composite surface is R* — R/y/2. 

For a given normal load on the solids, we note h the distance between the mean 
plane of the rough surface and the reference flat. The compression of a contacting 
asperity of height z > h is thus 5 = z — h (Greenwood & Williamson 1966). 
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(c) First loading of the multicontact interface 



We now calculate the interface response when the shear displacement x first 
increases from to x max . 

For each microcontact, the Hertz radius a and the normal load w depend on the 

compression 5 = z — h: 

a = ^/R* (z - h) (A3) 

E ^{z-hf 2 (A 4) 



AE* 3 _ 2 
W ~3lF a ~ 3(l-i/ 2 )" 



Plugging (A3) and (A 4) into equation (A 1), and using G = E/(2(l + u)), we 
obtain: 

2/3 



f 1 --)' 



1 



jj,\(z — h) 



where \ = s(2 - v)/[2(l - v)]. 

We note that the rhs term of equation (A 5) is > for: 

xs 

z-h> — 
/iA 



(A 5) 



(A 6) 



This means that microcontacts whose compression z — h satisfies condition (A 6) 
bear a tangential force / < \iw. Microcontacts such that z — h = xs/ (/iA) are totally 
sliding, i.e. f = fiw. 

In the following, we will assume that 'small' contacts such that z — h < xs/(fj,X) 
are also sliding and bear a tangential force equal to fiw]. 

Hence, for microcontacts such that (A 6) is verified, equation (A 5) leads to: 



2E 



f___ 

(i 3(1 - v 2 



nj*(z-k) 



3/2 



/tA(z — h) 



3/2' 



(A 7) 



The total tangential force on the system is obtained by integration over the 
height distribution: 



2E 



H 3(1-^) 



/ R*N, 



fro 



(z-hy 



3/2 



/tA / 



1 



' z/s dz 



+ J (z - hf /2 -e~ z ' s dz \ 



(A 8) 



where N is the total number of asperities, and ho is given by condition (A 6). 
The first integral corresponds to the contribution of microcontacts whose response 
is given by (A 7), while the second term corresponds to totally sliding contacts. 

Setting £ = [z — h)/s, and noting that the total normal load on the interface 
reads, according to Greenwood: 



r+oo 

W = (4/3)E*^N s 3/2 e- h/s / £ 

S v 'Jo 



3/2„-S 



(A 9) 



Wo 



f Actually, when the interface is sheared, the smallest microcontacts are destroyed and replaced 
by new ones, the contribution of which we cannot calculate. We will see further that the interfacial 
response is not significantly affected by this assumption. 
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we finally get the following expression for the macroscopic shear force: 

F = nW (l- (A 10) 

Inverting (A 10) yields: 

, = - M Aln(l--^) (All) 

In the limit of small tangential displacements, to lowest order in x/fiX <C 1, 
F = Wx/X. We thus find the expression of the interfacial shear stiffness n — W/X, 
where A = s(2 — v)/[2(l — v)] is the elastic length. For numerical purposes, we 
will however not estimate A from this expression but will rather make use of its 
experimentally measured value A ~ 0.25 /im. 

Finally, evaluating, from equation (A 8), the relative contribution of totally slid- 
ing microcontacts, we find that they contribute less than 10% to the calculated shear 
forcef while the tangential displacement stays lower than 60 nm (which corresponds 
to reduced shear force amplitudes 7 ac < 0.2, see e.g. figure 4). Their contribution 
is at most 20% for x ~ 100 nm. 



(d) Unloading from (F maxi x max ) 

Let us now study the case where the tangential displacement of the slider is de- 
creased, after having reached a maximum value x max , corresponding to a maximum 
shear force F max . Two families of microcontacts must then be considered: 



(i) Microcontacts such that (z — h) < x max s/{nX) 

At the end of the first loading, these microcontacts are totally sliding. For x = 
Xmax, the shear force on one of them is f ma x — /««, and its response when x is 
decreased is given by (Mindlin et al. 1953, Johnson 1985): 



where b is the inner radius of the corresponding slip zone: b = a[l/2 + f /(2/j.w)] 1 / 3 . 
We thus obtain: 

1 I f \ _ 1 ( X max ~ X\)S fAlTl 

2 2fiwJ 2/iX(z -h) ( ' 
As in §c, this equation yields the following condition: 

z-h> {Xmax ~ X ^ )S (A 14) 

2/xA v ; 

Microcontacts whose compression satisfies condition (A 14) respond according to 
equation (A 13), while those for which z — h < (x max — x\)s/(2fj,X) are assumed 
to be totally sliding and bear a tangential force / = —fiw. 

f This estimate is done using fi = 0.49, the value which leads to the best fit of xo(j ac ) (sec 
figure 9). A higher value of fi would lead to an even weaker contribution of sliding contacts. 
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For the contacts which are not fully sliding during unloading, equation (A 13) 
yields: 



2E 



f_ = 

jJL 3(1 - V 2 



/ R* I 2 



z — h 



-i 3/2 



(z-h) 



3/2 



(A 15) 



(ii) Microcontacts such that (z — h) > x max s / (jj,X) 

The response of these microcontacts, which were not totally sliding when x 
reached its maximum value, is given by (Mindlin et al. 1953, Johnson 1985): 



3(2 — v)jj,w 



2 1 



8Ga 

From equation (A 5) we get: 



fmax f 
2flW 



2/3 



1 



fn 



fJ,W 



2/3 



1 



1 - 



2/3 



ji\{z — h) 



(A 16) 



(A 17) 



which, once plugged into (A 16), leads to the following expression for the shear force 
on one of these microcontacts: 



2E 



/ 

fj, 3(1 - i/ 2 ) 



<R* { 2 



yUA 



z — h 

3/2 



2/iA 

(z-hM 



3/2 



(A 18) 



(iii) Total shear force 

The macroscopic shear force when the interface is unloaded is the sum of the 
contributions (A 15), (A 18), and of that resulting from fully sliding microcontacts: 



F 

- = W 
I 1 



2K- 



2^A 



3/2 



-u- 



3/2 



3/2 



+ 



J 27Ta 



2/iA 



3/2 



£3/2 



(A 19) 



which yields: 



= 2e 2 ^ A — e ma — l 



Setting 7 = F/W yields: 



(A 20) 



x\ = x max + 2/iA In 1 + 



2/z 



(A21) 
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with x max as given by (All): 

z max = -Mln(l-^p) (A 22) 

Finally, by symmetry, the relationship x(F) when the displacement, having 
reached a minimum value x m im is reversed up to x max , reads: 

r = .)■„,;„ - 2//A 111 ( 1 - 7 2 ^""" ) ( \ 2:-! 



with the minimum value of the displacement: 

X ml n = x max + 2fi\ In [1 I (A 24) 

With a macroscopic shear force of the form F(t) /W = -f(t) = + "f ac cos (cut), 
we obtain for the displacement response: 



x N = -/iAln ( 1 - I + 2^Aln 



7 ac (cosjejt) - 1) 
2pt 



(A 25) 



and 



X/ = -n\ha. ( 1 - 2^ ) +2/iAIn( 1- — I - 2^Aln 



_^ _ 7 ac (cosjeut) + 1) 
2^i 

(A 26) 



(e) In-phase and out-of-phase components 
The elastic and dissipative responses are given by: 



x 



x cos(u>t)dt 



x\ cos(ut)dt 



x / cos(uit)dt > (A 27) 



10 

v = — y a:\sm 



(ujt)dt + j 



x y sm(ojt)dt 



(A 28) 



We expand in power of j ac / (2/i) <C 1 the time-dependent logarithmic terms in 
equations (A 25) and (A 26), and finally get for the in-phase amplitude: 



xq = 2/j.X 



lac { lac X 5 / 7ac \ \ Q ( T ac 

2[i \2 f iJ 4:\2n) \2^ 



(A 29) 



The out-of-phase component can be calculated exactly and reads: 



^90 = 



4/xA 



1 _ 2E ) l n ( i _ 2^1 | - 2 



(A 30) 
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